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LONG-TERM ORBITAL EVOLUTION OF SHORT-PERIOD COMETS FOUND IN PROJECT ’’COSMO-DICE” 

Tsuko NAKAMURA (National Astronomical Observatory, Mitaka, 181 Tokyo) and 
Makoto YOSHIKAWA (Communications Research Laboratory, Kashima, Ibaraki, 314 Japan) 

ABSTRACT. Orbital evolutions of about 160 short-period (SP) comets are numerically 
integrated for 4400 years in the framework of a realistic dynamical model. By the 
round-trip error in closure test, a reliable time span of the integrated orbits is 
estimated for each comet. Majority of the SP comets with their Tisserand’s constant(J) 
between 2.8 and 3.1 are found to evolve within the past 1000-2000 years from the 
orbits whose perihelia are near the Jovian orbit to the orbits with perihelia of 1-2 
AU. This evolution is much more rapid than that expected from Monte Carlo simulations 
based on symmetric distribution of planetary perturbations, thus suggesting that 
asymmetry of perturbation distribution play an important role in cometary evolution. 
Several comets are shown to evolve from the near-Saturn orbits and then to be handed 
over under the control of Jupiter. We also find that a few comets were captured from 
long-period orbits (a = 75-126 AU) via only a few close encounters with Jupiter. It is 
confirmed that the captured SP comets of low-inclination with 2.8 < J < 3. 1 show more 
or less strong chaotic behavior. On the other hand, comets with longer orbital period 
and/or of high inclination reveal slow or quasi-per iodic orbital evolution. 

1. Introduction 

"Cosmo-DICE” is a project, named after Dynamical Investigation of Cometary Evolution 
and after the erratic nature of SP comet motions just like casting dice. By this 
project we intend to provide reliable results of long-term orbital evolution of SP 
comets in the framework of a dynamical model as realistic as possible both for comets 
and planets. Although this work may look like only a successor of Carusi et al. 
(1985a), there are some important differences between theirs and ours. First, the time 
span covered by our calculations is 8.5 times toward the past and 2.5 times toward the 
future longer than those by Carusi et al. (1985a), which are both 400 years. For this 
achievement, some specific devices described below were necessary in terms of choice 
of integrators and calculation accuracy. As a result of this, we could find 
substantial orbital evolution of many SP comets, whereas, in Carusi et al. (1985a), 
only very limited cases of orbital evolution can be seen. Second, our calculations are 
characterized by an elaborate error check of integration for each orbit. This enables 
us to have a reliable time interval of integration for each comet. 

Of course, considering the errors of the adopted orbital elements and possible non - 
gravitational effects which are not taken into account here, we cannot claim that 
every orbit in our calculations represents the true evolution for the 4400-year full 
time span. The main purpose of this project is, rather, to provide a standard database 
for statistical studies of long-term orbital evolution of SP comets, such as chaotic 
nature of thier orbits, capture into and ejection out of resonant motions, 
interrelation between asteroids and comets, and so on. 

In section 2 we describe briefly the dynamical model and method of integration. 
Section 3 discusses the reliable time span of our integration in relation to the error 
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growth rate. In section 4 are given the general characteristics of the orbital 
evolution of SP coaets found in this project. The details of orbital evolution for 
each coaet have been published in Nakaaura and Yoshikawa (1991). A FITS-foraatted 
aagnetic tape including all the results is also available on request fron the authors. 


2. Dynaaical Model and Method of Integration 

Orbital eleaents as the initial conditions for nuaerical integration are tak en f or 
the aost part froa the Catalogue of Coaetary Orbits, 6th edition (Narsden 1989) and 
the Minor Planet Circular. Planetary positions froa Mercury through Pluto, are not 
integrated simultaneously with coaets, but are taken froa the JPL nuaerical 
epheaerides DE102 (Newhall et al. 1983), which covers froa BC 1411 through AD 3002. 

In order to overcoae the orbital instability inherent in SP coaets, we adopted a 
quadruple-precision version of a variable-step extrapolation-type integrator (Bulirsch 
and Stoer 1966). The truncation paraaeter is set to 10" 22 at r = 1 All, which is a best 
coaproaise between CPU tiae and the accuracy requireaent froa orbital instability. 
Because of the variable-step control, regularization techniques are not used, 


3. Error Estiaate 

After several trials we found that the round-trip error in closure test (Carusi et 
al. 1985b) is best as an error checker of integration for our purposes. Naaely, 
integration is started froa the present epoch toward the past up to BC 1411, and then 
this final state is integrated back to the present. The difference of these two 
orbital paths at each tiae step is expected to be a good aeasure of the error in the 
nuaerical integration. This is because the integration toward the present should 
reproduce exactly the original initial state if there is no error in the integration. 

Fig. 1 shows soae exaaples in the tiae history of the absolute value of the round- 
trip error in true longitude. General trend of these curves is always of linearly 
growing nature on the seai-log scale, though the error growth rate can be different by 
several order of aagnitude froa coaet to coaet. The error g row th rat e is found to have 
intiaate relations with the orbital instability. This is reas ona ble because the 
Lyapunov characteristic exponent for a chaotic orbit should be positive (Lichtenberg 
and Lieberaan 1983) and thus the I inear error growth 1 ike Fig. 1 is expected. 

fe adopt here as a practical aeasure of orbital instabi 1 ity the elapsed tiae since 
BC 1411 in the backward integration within which the absolute error in true longitude 
attains one deg arc. This liait of tiae (T 1 i a) is calculated for each coaet and the 
correlation of Tlia wit h J value is plotted i n Fig. 2, fe see that the Bean Ilia for 
low-inclination SP coaets is about 2000 years, and the low-incli nation coaets with 
their J near 3.0 show chaotic aotions whereas the high-inclination ones behave aore 
regularly. The Tlia is also useful to estiaate the reliable tiae interval of the orbit 
calculation for a given error of adopted orbital eleaents. The error growth rate for 
the Bean Tlia of 2000 years is about 10 per century. The relative error of best 
orbital eleae nts deterained froa observations, on the other hand, is currently 10’ 7 or 
so. This iaplies that for an orbital longitude error of one deg arc the reliable tiae 
span of integration is about 700 years. 
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Fig. l: Soie examples of the tiie history of the absolute value of the round 

trip error in true longitude. The ordinate is in the scale of log(deg). 
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Fig. 2: Tisserand( Jacobi) const, vs. T-liiit 

4. General Trends of Orbital Evolution 

The following is soie overall characteristics of orbital evolution of the SP coiets 
found in this project, le have to keep in lind that these conclusions are lore or less 
inevitabl ly statistical and the orbital evolution of individual coiets should not 
always be taken at their face value. 

(1) Majority of the SP qoiets that show drastic change of perihelion distance (q) have 
the Tisserand’s constant (J) of 2. 7 — ^3. 1 : 

Since the solar radiation lakes coiets luiinous. the heliocentric distance or 
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equivalently q is the lost iiportant factor for observability Conditions. Therefore we 
concentrate our concern on the orbital evolution of q and classify the degree of q- 
change according to the value of J. In Table 1 group A is defined by the coiets whose 
J is larger than 2.9 and whose q undergoes the change bewteen qi= 1 — 2AU and q u =4~5AU 
during the 4400 years. Ve designate, on the other hand, as group a the coiets whose J 
is in the saie range but whose q shows the variation of no lore than 1~2 AU. 
Similarly, other groups are also defined as given in Table 1. This classification of 
course cannot be so strict, because considerable nuaber of calculated coiets showed 
the J-change of 0. 05 — 6.1. Soie coiets are classified with an apostrophed letter as 
critical cases of large q change, where qi=l~2 Al) but q„=~3.5 AU. 

When percentage fractions are calculated for each group in Table 1, a half nuiber of 
the apostrophed coiets are included in the upper- case group a nd another half in the 
lower-case group. It is reiarkable to see that 85% of the coiets with 2.9 J < 3.1 
belongs to group A. This is partially explained by an approximate equation U 2 = 3 - J, 
where J is an expression of Tisserand’s invariant and U is the total encounter 
velocity with Jupiter under the assuiption of two-body interaction (*0*p i k 1963); for 
group A, U becoaes nearly zero, so that large perturbations can be induced. 

(2) The nuiber of the large q-change coiets whose q, seii-iajor axis, and aphelion 
distance decrease globally in their orbital evolution is iuch lore than that of the 
coiets that evolve in the opposite sense: 

Along the 85 capital-letter coiets, 59 (69%) showed the decreasing evolution, 9(11%) 
the increasing one, and 17 (20%) no specific trend. The decreasing evolution takes 
place lostly in the past, while the increasing one in the future. Therefore, in order 
to lake unbiased comparison between the decreasing and increasing evolutions, the tiie 
span ratio for the future and the past (3400/1000) lust be corredted to the nuiber of 
the increasing evolution coiets. However, the corrected nuiber of the increasing 
evolution coiets is 3.4 x 9 = 31, which is still only a half of the dereasing 
evolution coiets. This leans that capture process is lore favoured than ejection in 
the evolution of SP coiets. This finding is consistent with the reported asyntery in 
the distribution of energy perturbation (Nakaiura 1981; Nanara and Valsecchi 1991) 
which prefers the paths toward shoter-per iod orbits. Interestingly, on the contrary, 
Everhart(1969) and Carusi et al . ( 1990) report the opposite tendency. This discrepancy 
will thus be a subject of f urther investigations. 

(3) Most of the orbits^ before capture into Jupiter family have q ~ 5 AU: 

All of the 59 captured SP coiets lentioned in the previous paragraph were originated 
froi the orbit s of q ~ 5 AU. These coie ts correspond to the ones in the "capture 
region” by Jupiter, first confined numerically by Everhart (1972). The orbits before 
capture are also characterized by the very flat nature of q. Eight out of the 59 
captured coiets are shown to evolve froi the orbits near Saturn and then to be relayed 
under the control of Jupiter; this is the first aanifestation of the lultiple-stage 
capture aechanisi proposed by Everhart (1977). Soie of the exaiples are given in Fig. 
3. P/He 1 in-Ro aan- Alu 2, P/Honda-Mrkos-Pajdusakova, and P/Barnard 3 are found to be 
captured froi the iong -period orbits of a =75~ 125 AU into SP orbits through a few 
close encounters with Jupiter (Fig. 4). . - • 

(4) The evolutionary tiie scale for capture is much shorter than that predicted by 
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Monte Carlo simulations of simplified dynamical models* 

We found many cases of the orbital evolution in which the orbital period (P) and q 
jumped from P=12~20 yr and q— 5 AU into P=6— 7 yr (typical of SP comets) and 0=1—2 
AU respectively via only one or two strong encounters with Jupiter. In analytical 
theory and Monte Carlo simulations based on symmetric distribution of perturbation, 
drastic orbital changes in such a short-time scale are very rare (e.g., Yabushita 
1980; Nakamura 1981). This suggests that the asymmetric structure in the tail of 
perturbation distribution play an important role for the capture into SP comets. 

Table 2 is a distribution of how many years ago an abrupt capture took place for the 
59 comets analyzed above. It is seen that 
nearly one third of the comets were cap- 
tured very recently (within the past 200 
years) and they have been discovered soon 
after their capture. This may suggest that 
overlooked comets in observations be fewer 
than used to be anticipated as observational 
selection. 


References 

Bulirsch R. and Stoer J . (1966) Numerische Mathematik, 8, 1-13. 

Carusi A., Kresak L., Perozzi E., and Valsecchi G.B. (1985a) Long-term orbital 
evolution of sh ort-period comets, Adams Hilger (Bristol). 

Carusi A., Kresak L., Perozzi E. , and Valsecchi G.B. (1985b) in Dynamics of Comets: 
Their origin and evolution, eds. by A. Carusi and G.B. Valsecchi, p203— 2 14, Reidel 
(Dordrecht). 

Carusi A., Valsecchi G.B., and Greenberg R,(1990) Celestial Mechanics and Dynamical 
Astronomy, 42, 111-131. 

Everhart E.(1969) Astron. J., 74, 735-750. 

Everhart E. ( 1972) Astrophys. Letters, JO, 131-135. 

Everhart E . ( 1977 ) in Comets Asteroids Meteorites, ed. by A.H. Delsemme, p99- 1 04 , The 
University of Toledo. 

Lichtenberg A.J. and Lieberman M. A. (1983) Regular and stochastic motion, Springer- 
Verlag (New York). 

Manara A. and Valsecchi G.B. (1991) Abstract of this symposium. 

Marsden B.G.(1989) Catalogue of cometary orbits, 6th edition, Central Bureau for 
Astronomical Telegrams. 

Nakamura T. (1981) Icarus, 45, 529-544. 

Nakamura T. and Yoshikawa M . (199 1 ) Pulb. Natl. Astron. Obs. Japan, 2» 293-383. 

Newhall X.X., Standish Jr.E.M., Williams J.G. (1983) Atron. Astrophys., 84. 1910-1917. 
Op ik E. J . ( 1963) Advances Astron. Astrophys., 2. 219-262. 

Yabushita S . ( 1 980) Mon. Not. Roy. astr. Soc., 190, 71-73. 


Table 2 - Distribution of capture time 


Interval (yr) 

Number 

0 

- -200 

1? 

-200 

- -500 

9 

-500 

- -1000 

5 

-1000 

— -1500 

n 

-1500 

- -2000 

9 

-2000 

— -3400 

8 



